Magnetic properties of LaOi-^F^FeAs 
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Using state-of-the-art first-principles calculations we have elucidated the complex magnetic and 
structural dependence of LaOFeAs upon doping. Our key findings are that (i) doping results in an 
orthorhombic ground state and (ii) there is a commensurate to incommensurate transition in the 
magnetic structure between x = 0.025 and x = 0.04. Our calculations further imply that in this 
system magnetic order persists up to the onset of superconductivity at the critical doping of x = 
0.05. Finally, our investigations of the undoped parent compound reveal an unusually pronounced 
dependence of the magnetic moment on details of the exchange-correlation (xc) functional used in 
the calculation. However, for all choices of xc functional an orthorhombic structure is found. 

PACS numbers: 74.25. Jb,67.30.hj,75.30.Fv,75.25.tz, 74.25. Kc 



The discovery of high-Tc superconductivity in the iron 
oxypnictide LaOi-^F^FeAs is attracting a lot of atten- 
tion, mostly because this class of materials is one of the 
first set of materials to have a high T c like Cuprates [HQ, 
with a striking difference: in oxypnictides the itinerant 
magnetic FeAs layer plays a strong role in superconduc- 
tivity, while in the cuprates this role is adopted by highly 
localized moments in the CuO plane. This of course puts 
a spotlight on the involvement of magnetic fluctuations 
0,0, SB] in the onset of the superconducting state, mak- 
ing iron oxypnictide a playground for understanding the 
mechanism and role of magnetism in high T c supercon- 
ductivity. 

Understanding the magnetic ground-state of 
LaOi-^F^FeAs is the first crucial step towards un- 
covering the mystery of the role of magnetism in the 
onset of superconductivity. Determining magnetic 
ground states using first principles density functional 
theory (DFT) calculations is a routine procedure. 
However, in the present case this has proved to be a very 
difficult task. All experiments, whether probing long 
range order via neutron scattering @, H i, E3] or local 
probe studies via the Mossbauer effect 9( , conclusively 
report an itinerant small (0.25-0.36 /ib) moment per Fe 

atom @,H,[l,[lQj]. 0n tne otner hand; DFT calculations 
(performed using both experimental and optimised 
atomic positions) have resulted in moments ranging 
from 0A7fiB to more than 2.0/^s. This dramatic spread 
of results is presumably indicative of both a sensitive 
dependence of the magnetism upon calculational details, 
as well as a strong dependence on crystal structure. 

In fact, one of the most complex and interesting as- 
pects of magnetism in LaOi^xF^FeAs is its interplay 
with structural properties. At a temperature of 136 K the 



undoped parent compound undergoes a phase transition 
to a stripe antiferromagnetic (AFM) order which, addi- 
tionally, is associated with a structural distortion from 
tetragonal to orthorhombic fH lil f]~oj] . The impact of dop- 
ing upon structural and magnetic properties is only be- 
ginning to be understood fa. : experiments differ on 
whether doping (greater than 5%) suppresses the struc- 
tural phase transition. Furthermore, a recent experiment 
finds evidence of incommensurate order upon doping ll| , 
indicating a rich interaction between structural and mag- 
netic order as the parent compound is doped. Highly 
accurate ab-initio calculations can play a crucial role in 
elucidating this complex behaviour, and in the present 
work we shall address in detail the magnetic and struc- 
tural properties of both doped and undoped LaOFeAs. 

In order to keep the numerical analysis as accurate 
as possible, in the present work all calculations are per- 
formed using the state-of-the-art full-potential linearized 
augmented plane wave (FPLAPW) method [l2|, imple- 
mented within the Elk code fl3| . We have taken great 
care that all relevant calculational parameters are con- 
verged; in particular we use a k mesh of 12 x 12 x 8 shifted 
by [0.5,0.5,0.5], and 188 states per k-point which ensures 
convergence of the second variational stepJUJ. We found 
that using less well converged values of these parameters 
has a dramatic impact on the magnetic moment (such 
computational details can be seen in the supplementary 
materail). 

We first determine the magnetic ground state of 
the undistorted (tetragonal phase) parent compound 
LaOFeAs. That magnetism in this system is remarkably 
sensitive to calculational details is clear from the spread 
of moments obtained using the FP-LAPW method and 
experimental crystal structure, 0.87-2.20/j,b 3, 15, 0, 
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FIG. 1: (color online) Top-left hand panel shows the magnetic 
moment per Fe atom (in \ib) and lower- left hand panel the 
total energy (in meV) per formula unit. All quantities are 
plotted as a function of spin-spiral (/-vector. The right hand 
panel displays the magnetic moment per Fe atom (in ^b) as a 
function of position of the As atom in the stripe AFM phase. 
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In this regard it has been noted that both 
Fermiology [l6| as well as magnetic moment [l9[ are 
strongly dependent on the z coordinate of the As atom. 
In the right hand panel of Fig. [T]we show this dependence 
of the magnetic moment on the z coordinate of the As 
atom for the so-called stripe AFM phase EH, Ei (the 
proposed ground state structure of this system), for two 
common parameterizations [20I [2l| of the local spin den- 
sity approximation (LSDA). Clearly, this choice leads to 
a dramatic difference in the results. For instance, at the 
value of zas obtained by optimizing with the generalised 
gradient approximation (GGA) (z^ A = 0.638 a. u.), the 
difference in the moment per Fe atom calculated using 
two different parameterizations for LSDA is 0.44/xb- It 
is further clear that only for a theoretically optimised za s 
(GGA or LSDA) can a moment close to the experimen- 
tal value be obtained. Taken together with the work of 
Mazin et ai, where it was shown [l6| that only for z^ A 
a Fermiology consistent with experiment could be found, 
these results indicate that a theoretical value of this pa- 
rameter should be used. However, given this unusually 
strong dependence of magnetism upon the choice of ex- 
change correlation functional, we adopt the strategy of 
calculating structural properties with the GGA, but then 
calculate ma gne tic properties via both the von Barth- 
Hedin (vBH)0 and Perdew-Wang (PW)@ parama- 
terizations of the LSDA. 

Given this, we now explore the phase space of pos- 
sible magnetic structures consistent with experiments. 
To date, the commensurate collinear stripe and checker- 
board structures have been investigated, with the stripe 
phase substantially lower in energy EU 22 1. Two cru- 



stripe phase remain the minimal energy structure? and 
(ii) what is the dependence of the moment on the under- 
lying spin structure? 

We here calculate spin spirals with q around the 
[1/2,0,0] (X) and [1/2,1/2,0] (M) special points in the 
Brillouin zone. Two special cases of these are the com- 
mensurate stripe and checker-board phases, generated 
with wave vectors M and X, respectively. Around the 
M -point, an alternative magnetic order is a spin spiral 
with a 90° phase difference between the two Fe sites, 
equivalent to a 2k structure at the M-point. Due to 
the inherent frustration of ordering at the M-point, this 
structure is also a candidate for the ground state [23| . 
In the upper-left panel of Fig[T] are shown the moment 
per Fe atom for spin spiral vectors with q near the M 
point for lk type spin spiral. Remarkably, we find that 
for the vBH parameterization a magnetic solution exists 
only around the M point and even a small change in spi- 
ral q from [0.50,0.50,0.00] to [0.62,0.62,0.00] results in 
a vanishingly small value of the moment (similar results 
were obtained using LMTO method in Ref. 24). The 
PW parameterization also results in a moment which falls 
stongly on moving away from the M point, but the value 
is non-zero over the whole BZ with a minimal value of 
0.20fiB at the T point. 

Turning to the energetics of these spirals (lower-left 
panel FigfT]) , we find the lk type spiral to be always lower 
in energy than the 2k type spiral (for clarity we show only 
the 2k result for vBH-LSDA). Interestingly, the global en- 
ergy minima is not at the M point but instead at a spiral 
of q = [0.52, 0.52, 0.00], a results which is independent of 
the functional chosen. The ground state of undistorted 
LaOFeAs is thus an incommensurate spin spiral. With 
the vBH LSDA this structure lowers the total energy by 
0.5 meV per formula unit as compared to the previously 
supposed stripe ground state TH, 1|| and gives a moment 
of 0.46/is- This energy difference and the moment are 
increased to 2 meV per formula unit and 0.92/1^ respec- 
tively when the PW LSDA is used for determining both 
the zas an d the magnetic structural energies. 

At a temperature of 150 K LaOFeAs undergoes a phase 
transition from the tetragonal to orthorhombic structure, 
closely associated with the onset of magnetism[3, 0, EJ 
13, EE S3, Si. Recently, this transition was explained as 



cial questions which then may be asked are: (i) when 
the constraint of commensuration is removed does the 



the system lowering magnetic frustration, rationalised by 
arguments based on a J\ — J2 Heisenberg model. Here we 
shall take a different aproach and explain the mechanism 
of this distortion without recourse to the J\ — J2 model, 
which is unlikely to be appropriate for a weak itinerant 
magnet such as LaOFeAs. 

In Fig. [2] (upper panel) is displayed the total energy 
of the non-magnetic state of LaOFeAs as a function of 
distortion angle, for both the vBH and PW parameteri- 
sations; clearly, a special feature is the very flat minima 
at 90°. Considering first the commensurate stripe phase 
we find that the moment (middle panel) increases with 
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FIG. 2: (color online) Shown are the total energies for the non- 
magnetic system (top panel) and the magnetic stripe phase 
(bottom panel), in meV per formula unit. In the middle panel 
is displayed the corresponding moment of the stripe phase (in 
\ib per Fe atom). All quantities are plotted as a function 
of distortion angle; 90° corresponds to undistorted unit cell. 
Insets in the bottom panel illustrate the distorted magnetic 
structures. 



increasing distortion angle. The gain in magnetic energy 
is then sufficient to shift the minimum of the total energy 
from 90° to 90.30° (lower panel) for the vBH and 90.50 
for the PW LSDA, in good agreement with the mea- 
sured value of 90.27°. Thus the vBH LSDA gives both 
the ground state moment and distortion in close agree- 
ment with experiment. On the other hand, the larger 
magnetic moment produced by the PW parameterisation 
simply results in a slightly larger distortion; reasuringly 
therefore, the qualitative behaviour is independent of the 
particular choice of functional. 

Turning now to possible non-collinearity of spins in this 
distorted phase of LaOFeAs we find that, in contrast to 
the undistorted crystal, the commensurate stripe phase 
is the minimal energy structure in the distorted case (for 
both parameterizations of LSDA) . We should stress that 
the analysis of the structural distortion pesented here 
relies solely on the behaviour of the global moment with 
distortion and the concomitant increase in magnetisation 
energy, and is thus appropriate for the itinerant nature of 
LaOFeAs. In Ref.llathe interesting suggestion was made 
that the cause of the distortion could be a release of the 
inherent magnetic frustration of the FeAs layers. Such a 
relaxation of frustration would, generally, be expected to 
lead to an increase in moment, and so this idea is also 
compatable with the picture presented here. 

The electron doped material LaOi-^F^FeAs becomes 

0,1" 



superconducting [Ij, |8|, |ll|, |27j at a critical doping of 
x = 0.05. A crucial question is whether magnetism 
persists up to the superconducting transition, i.e., is a 
competing ground state, or is lost before the onset of 
superconductivity. Experimentally, it is established that 
beyond x — 0.075 no magnetic order persists, however if 
magnetic order is entirely lost before the superconduct- 
ing transition is still a point of discussion [§[ 11, 28|. It 
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FIG. 3: (color online) Spin magnetic moment per Fe atom (in 
[ib) as a function of electron doping. The experimental data 
is taken from Ref. [Til . 



is also evident from recent experiments [ll| that there 
exists a complex structural and magnetic behaviour with 
doping, the full nature of which has yet to be clarified. In 
fact, on the question of whether the distortion observed 
at x = persists beyond the superconducting transition 
at x = 0.05 , the current experimental data are contra- 
dictory; Rcfs. 



11 and 



J find that distortion persist up 
to x = 0.08, while in Refs. H, no distortion is observed 
beyond the onset of superconductivity at x = 0.05. The 
magnetic state of of LaOi-^F^FeAs is also uncertain, 
with one experiment ll| finding evidence of an incom- 
mensurate structure for x > 0. 

In order to clarify this situation we have determined 
the ground state for several doping concentrations, by 
minimising over both the distortion angle 7, see Fig. 
|4j and spiral vector q (for lk type spin configuration). 
To simulate doping a small amount of charge is added 
to the unit cell along with a compensating background, 
which ensures charge neutrality. Since the concentration 
of electrons is very small this approximation to the dop- 
ing should be good. Indeed, the accuracy of this approx- 
imation has recently been demonstrated [2^]. 

We first consider the overall behaviour of the magnetic 
moment with doping. This is shown in Fig. [3] for both 
the vBH and PW parameterisations of the LSDA. It is 
immediately apparent that the choice of functional sim- 
ply leads to a scaling of the doping curve. However, only 
for the vBH LSDA do we find a behaviour that results in 
quantatitive agreement with available experimental data 
for x < 0.05, i.e., less than the critical doping. How- 
ever, a marked divergence between experiment and the- 
ory occurs after this point; while the theoretical data 
show a slowly vanishing tail the experimental data dis- 
play a sudden decrease. Our calculations, therefore, im- 
ply that magnetic and superconducting order are com- 
peting ground states in the sense that magnetism does 
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not die before the onset of superconductivity. The sharp 
decrease in moment at x = 0.05 may then be brought 
about by the onset of superconductivity. The question of 
whether magnetism and superconductivity then coexist 
for a small doping range, or if the onset of superconduc- 
tivity destroys entirely the magentic order, we cannot 
answer. On the other hand, one should note that the 
small moment itinerant magnetism of this system im- 
plies a strong role for spin fluctuations [3, @, Ho] which is 
not correctly treated by the LSDA functional used in the 
present calculations. These may act to damp the slowly 
vanishing tail of the moment vs. doping seen in Fig. [3] 




mensurate and stripe structures (0.4 meV per formula 
unit) suggests that this may only be seen at the low tem- 
perature (8 K) experiments performed by Huang et al. 
Furthermore, in experiments a small increase in the pa- 
rameter 2^ pt is seen with doping, an effect we have not 
included in our calculations. As increasing this param- 
eter results in an increase of moment, the likely effect 
of including this would be to shift the onset of moment 
lowering mechanisms (7 < 90° distortion and incommen- 
surate order) to somewhat higher values of doping. 

To summarize by the means of accurate ab-initio calcu- 
lations we have given an explanation for the phase transi- 
tion in LaOFeAs based on an itinerant magnetic picture 
in terms of an increase in spin polarisation brought about 
by crystal distortion. Furthermore, we have elucidated 
the impact of doping on the ground state; we find both 
that the doped material distorts, but with 8 = 89.85, and 
a stripe phase to incommensurate spin spiral transition 
takes place between x — 0.025 and x = 0.04. Most im- 
portantly, our calculations indicate that in this system 
magnetic order persists up to the onset of superconduc- 
tivity. 

Supplementary material 
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FIG. 4: (color online) Left panel shows the total energy (in 
meV per formula unit) as a function of distortion angle, for 
various dopings. Right panels show total energy (in meV per 
formula unit) as a function of spin-spiral q-vector for distorted 
(upper panel) and undistorted (lower panel) doped LaOFeAs. 

We now turn to the important question of the impact 
of doping upon the crystal and magnetic ground state. 
Given that the vBH parameterised LSDA has shown good 
agreement with experiment for the magnetic moment, 
distortion and doping curve, we shall use this functional 
to explore this question. We first fix the magnetic struc- 
ture to that of the commensurate stripe phase and, for 
each value of x, minimise the distortion angle. Remark- 
ably, we find for all x > a crystal distortion opposite 
(7 < 90°) to that of the undoped parent compound (see 
Fig. Hleft panel). Furthermore, subsequent minimisation 
of the spin spiral q reveals a stripe phase to incommen- 
surate spin spiral transition (q = [0.54,0.54,0]) between 
x = 0.025 and x = 0.04 (top right panel, Fig. 3). This 
opposite distortion upon doping can be understood as a 
mechanism the system adopts to lower the moment. The 
monotonic increase of the moment with 7, displayed in 
Fig. [2] for x — 0, is also found for all doping concentra- 
tions, and hence a crystal distortion of 7 < 90° lowers 
the moment. 

This finding of a commensurate stripe to incommen- 
surate transition is in agreement with the experimental 
observation [llj of Huang et al.. We should note, how- 
ever, that the small energy difference between the incom- 
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FIG. 5: (color online) Left panel shows the change in moment 
per Fe atom (in /is) as a function of number of k-points in 
the IBZ, for these calculations number of states per k-point 
is fixed to 148. Right panel shows the change in moment per 
Fe atom (in ub) as a function of the number of states per k- 
point, for these calculations number of k-points is kept fixed 
at 312 (unshifted mesh). 

In the present work all calculations are performed using 
the state-of-the-art full-potential linearized augmented 
plane wave (FPLAPW) method. To obtain the Pauli 
spinor states, the Hamiltonian containing only the scalar 
potential is diagonalizcd in the LAPW basis: this is the 
first- variational step. The scalar states thus obtained are 
then used as a basis to set up a second- variational Hamil- 
tonian with spinor degrees of freedom. This is more effi- 
cient than simply using spinor LAPW functions, but care 
must be taken to ensure that there is a sufficient num- 
ber of first- variational eigenstates for convergence of the 
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second- variational problem. In Fig. [5] are shown the re- 
sults for magnetic moment per Fe atom, for undistorted 
undoped LaOFeAs as a function of the number of these 
first- variational eigenstates per k-point. The results are 
for the stripe anti-ferromagnetic phase. It is clear that 
ground state magnetic moment for this material is greatly 
sensitive to both the number of states per k-point as well 
as the number of k-points in the irreducible Brillouin- 
zone (IBZ). 




FIG. 6: (color online)Top panel shows the magnetic moment 
per Fe atom (in /is), calculated using various methods, with 
As atoms at experimental position. Bottom panel shows the 
same but for As atom at theoretically optimized position, this 
position is obtained using GGA functional and by treating 
LaOFeAs non-magnetically. 

The magnetism in this system is remarkably sensitive 
to calculational details is clear from Fig. [6l In this fig- 
ure are shown the results of magnetic moments obtained 
using various approximations to the exchange-correlation 
functionals implemented within different full-potential as 
well as pseudo potential codes. The moments obtained 
using different methods and experimental As position, 
range from 0.87-2.20/xs and the moments obtained us- 
ing the theoretically optimized (GGA) As position range 



from 0.46-1.80/ib. On the other hand both Neutron 
diffraction and Mossbauer experiments report a value be- 
tween 0.25-0.36/i_B 
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